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We give an explicit expression for the geometric measure of entanglement for three qubit states that are linear 
combinations of four orthogonal product states. It turns out that the geometric measure for these states has 
three different expressions depending on the range of definition in parameter space. Each expression of the 
measure has its own geometrically meaningful interpretation. Such an interpretation allows oneself to take one 
QO . step toward a complete understanding for the general properties of the entanglement measure. The states that lie 

on joint surfaces separating different ranges of definition, designated as shared states, seem to have particularly 
<T^> interesting features. The properties of the shared states are fully discussed. 
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I. INTRODUCTION 



Entanglement is the most intriguing feature of quantum mechanics and a key resource in quantum information science. One 
, of the main goals in these theories is to develop a comprehensive theory of multipartite entanglement. Various entanglement 
measures have been invented to quantify the multi-particle entanglement [2, HlBl but none of them were able to suggest 
Qh| a method for calculating a measure of multipartite systems. This mathematical difficulty is the main obstacle to elaborate a theory 
■ of multi-particle entanglement. 

In this paper, we present the first calculation of the geometric measure of entanglement 10, [H, [H] for three qubit states which 
are expressed as linear combinations of four given orthogonal product states. Any pure three qubit state can be written in terms 
Q-i of five preassigned orthogonal product states floll via Schmidt decomposition. Thus the states discussed here are more general 
states compared to the well-known GHZ ifTTll and Wld states. 

The reason for using the geometric measure of entanglement is that it is suitable for any partite system regardless of its 
dimensions. However, analytical computation for generic states still remains as a great challenge. The measure depends on 
entanglement eigenvalue A^ ax and can be derived from the formula E g (ip) = 1 — A^ lax . For pure states, the entanglement 
eigenvalue is equal to the maximal overlap of a given state with any complete product state. This measure has the following 
remarkable properties: 

i) it has an operational treatment. The same overlap A^ ax defines Groverian measure of entanglement lfl3l [Till which has 
been introduced later in operational terms. In other words, it quantifies how well a given state serves as an input state to Grover's 
search algorithm 01511 . From this view, Groverian measure can be regarded as an operational treatment of the geometric measure. 

OO ' j — L 

li) it has identified irregularity in channel capacity additivity B16I1 . Using this measure, one can show that a family of quantities, 

which were thought to be additive in an earlier papers, actually are not. For example, it is natural to conjecture that preparing 

two pairs of entangled particles should give us twice the entanglement of one pair and, similarly, using a channel twice doubles 

| its capacity. However, this conjecture claiming additivity has proved to be wrong in some cases. 

iii) it has useful connections to other entanglement measures and gives rise to a lower bound on the relative entropy of 
entanglement |[T7ll and generalized robustness 018ll . For certain pure states the first lower bound is saturated and thus their 
relative entropy of entanglement can be deduced from their geometric measure of entanglement. The second lower bound to 
generalized robustness can be express in terms of A^ ax directly. 

Owing to these features, the geometric measure can play an important role in the investigation of different problems related to 
entanglement. For example, the entanglement of two distinct multipartite bound entangled states can be determined analytically 
in terms of a geometric measure of entanglement 01911 . Recently, the same measure has been used to understand the physical 
implication of Zamolodchikov's c-theorem [20] more deeply. It is an important application regarding the quantum information 
techniques in the effect of renormalization group in field theories 12111 . Thus it is natural that geometric measure of entanglement 
is an object of intense interest and in some recent works revised 12211 and generalized [23] versions of the geometric measure 
were presented. 

The progress made to date allows oneself to calculate the geometric measure of entanglement for pure three qubit systems 
rt24ll . The basic idea is to use (n — l)-qubit mixed states to calculate the geometric measure of n-qubit pure states. In the case 
of three qubits this idea converts the task effectively into the maximization of the two-qubit mixed state over product states 
and yields linear eigenvalue equations [25]. The solution of these linear eigenvalue equations reduces to the root finding for 
algebraic equations of degree six. However, three-qubit states containing symmetries allow complete analytical solutions and 
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explicit expressions as the symmetry reduces the equations of degree six to the quadratic equations. Analytic expressions derived 
in this way are unique and the presented effective method can be applied for extended quantum systems. Our aim is to derive 
analytic expressions for a wider class of three qubit systems and in this sense this work is the continuation of Ref. lE5ll . 

We consider most general W-type three qubit states that allow to derive analytic expressions for entanglement eigenvalue. 
These states can be expressed as linear combinations of four given orthogonal product states. If any of coefficients in this 
expansion vanishes, then one obtains the states analyzed in B25I1 . Notice that arbitrary linear combinations of five product 
states [10] give a couple of algebraic equations of degree six. Hence Evariste Galois's theorem does not allow to get analytic 
expressions for these states except some particular cases. 

We derive analytic expressions for an entanglement eigenvalue. Each expression has its own applicable domain depending 
on state parameters and these applicable domains are split up by separating surfaces. Thus the geometric measure distinguishes 
different types of states depending on the corresponding applicable domain. States that lie on separating surfaces are shared by 
two types of states and acquire new features. 

In Section II we derive stationarity equations and their solutions. In Section III we specify three qubit states under consid- 
eration and find relevant quantities. In Section IV we calculate entanglement eigenvalues and present explicit expressions. In 
Section V we separate the validity domains of the derived expressions. In Section VI we discuss shared states. In section VII we 
make concluding remarks. 



II. STATIONARITY EQUATIONS 

In this section we briefly review the derivation of the stationarity equations and their general solutions 1I25I1 . Denote by p ABC 
the density matrix of the three-qubit pure state and define the entanglement eigenvalue A, 2 nax ]7|] 

A^ax = m t x , tr (t jABO Q x ® e 2 ® g 3 ) , (i) 

gl gZ gi 

where the maximization runs over all normalized complete product states. Theorem 1 of Ref. [24] states that the maximization of 
a pure state over a single qubit state can be completely derived by using a particle traced over density matrix. Hence the theorem 
allows us to re-express the entanglement eigenvalue by reduced density matrix p AB of qubits A and B 



A4 ax = maxtr(p A V®^). (2) 

Now we introduce four Bloch vectors: 

1) ta for the reduced density matrix p A of the qubit A, 

2) r B for the reduced density matrix p B of the qubit B, 

3) u for the single qubit state g 1 , 

4) v for the single qubit state g 2 . 

Then the expression for entanglement eigenvalue (f2]l takes the form 

A roax = 7 , ma , x (1 + u ■ r A + v ■ r B + gijUiVj) , (3) 

4 U 2 — V 2 — 1 

where(summation on repeated indices i and j is understood) 

9ij =tT(p AB a i ®<j j ) (4) 
and <7,'s are Pauli matrices. The closest product state satisfies the stationarity conditions 

r A + gv = X 1 u, r B +g T u = \ 2 v, (5) 
where Lagrange multipliers Ai and A2 enforce the unit Bloch vectors u and v. The solutions of Eq.© are 

u = (AiA 2 l - gg T ) 1 (A 2 r\4 + gr B ) , v = (AiA 2 l - g T g) 1 (\\r B + g T r A ) . (6) 
Unknown Lagrange multipliers are defined by equations 
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u 2 = 1, v 2 = 1. (7) 

In general, Eq.© gives algebraic equations of degree six. The reason for this is that stationarity equations define all extremes 
of the reduced density matrix p AB over product states, regardless of them being global or local. And the degree of the algebraic 
equations is the number of possible extremes. 

Eq.© contains valuable information. It provides solid bases for a new numerical approach. This can be compared with the 
numerical calculations based on other technique l26ll . 

III. THREE QUBIT STATE 

We consider W-type state 

\-iJj) = a|100) + 6|010) + c|001> + djlll), (8) 

where free parameters a, b, c, d satisfy the normalization condition a 2 + b 2 + c 2 + d 2 = 1. Without loss of generality we consider 
only the case of positive parameters a, b, c, d. At first sight, it is not obvious whether the state allows analytic solutions or not. 
However, it does and our first task is to confirm the existence of the analytic solutions. 

In fact, entanglement of the state Eq.® is invariant under the permutations of four parameters a, b, c, d. The invariance under 
the permutations of three parameters a, 6, c is the consequence of the invariance under the permutations of qubits A,B,C. Now 
we make a local unitary(LU) transformation that relabels the bases of qubits B and C, i.e. Ob <-* 1b, Oc *-* 1(7 > an d does 
not change the basis of qubit A. This LU-transformation interchanges the coefficients as follows: a <-> d, b <-> c. Since any 
entanglement measure must be invariant under LU-transformations and the permutation b <-> c, it must be also invariant under 
the permutation a <-> d. In view of this symmetry, any entanglement measure must be invariant under the permutations of 
all the state parameters a, b, c, d. Owing to this symmetry, the state allows to derive analytic expressions for the entanglement 
eigenvalues. The necessary condition is |25| 

det(AiA 2 l -gg T ) = 0. (9) 

Indeed, if the condition (0 is fulfilled, then the expressions © for the general solutions are not applicable and Eq.© admits 
further simplification. 

Denote by i, j, k unit vectors along axes x, y, z respectively. Straightforward calculation yields 

(2lj \ 

r A = nk, r B =r 2 k, g=\0 2^ , (10) 

\0 -r 3 J 

where 

ri =b 2 +c 2 -a 2 ~d 2 , r 2 = a 2 +c 2 -b 2 ~d 2 , r 3 = a 2 + b 2 - c 2 - d 2 , (11) 
U> = ab + dc, fi — ab — dc. 

Vectors u and v can be written as linear combinations 

u = U{i + Ujj + Mfcfe, v = Vii + Vjj + Vkk (12) 

of vectors i, j, k. The substitution of the Eq.dT2T> into Eq.© gives a couple of equations in each direction. The result is a system 
of six linear equations 

2uiVi = X\Ui, 2u>Ui — \2Vi, (13a) 



2jjLVj — XiUj, 2/iUj = \2Vj, (13b) 
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ri - r 3 v k = AiUfe, r 2 - r 3 u k = A 2 u fc . (13c) 

Above equations impose two conditions 

(AiA 2 - 4uj 2 ) Ui Vi = 0, (14a) 

(A X A 2 - ifi^UjVj = 0. (14b) 

From these equations it can be deduced that the condition (O is valid and the system of equations © and (|7]i is solvable. Note 
that as a consequences of Eq.(fT3l x and/or y components of vectors u and v vanish simultaneously. Hence, conditions ( fT4b are 
satisfied in following three cases: 



• vectors u and v lie in xz plane 



vectors u and v lie in yz plane 



• vectors u and v are aligned with axis z 



AiA 2 -4cj 2 = 0, ii,-Uj=0, (15) 



AiA 2 - 4/i 2 = 0, = 0, (16) 



UiVi = UjVj = 0. (17) 

These cases are examined individually in next section. 

IV. EXPLICIT EXPRESSIONS 

In this section we analyze all three cases and derive explicit expressions for entanglement eigenvalue. Each expression has 
its own range of definition in which they are deemed applicable. Three ranges of definition cover the four dimensional sphere 
given by normalization condition. It is necessary to separate the validity domains and to make clear which of expressions should 
be applied for a given state. It turns out that the separation of domains requires solving inequalities that contain polynomials of 
degree six. This is a nontrivial task and we investigate it in the next section. 

A. Circumradius of Convex Quadrangle 
Let us consider the first case. Our main task is to find Lagrange multipliers Ai and A 2 . From equations d 1 3cb and ( fTBT ) we have 

A 2 ri-r 2 r 3 A x r 2 - rxr 3 

Auj z — 73 4lu z — r| 

In its turn Eq.( 113ab gives 

Xml = X 2 vl (19) 

Eq.(|7]l allows the substitution of expressions ( TT~8b into Eq.dT9b. Then we can obtain the second equation for Lagrange multi- 
pliers 

Ai (4^ 2 + r\ - rl) = A 2 (4^ 2 + r\ - r 2 ) . (20) 

This equation has a simple form owing to condition (O. Thus we can factorize the equation of degree six into the quadratic 
equations. Equations d20l ) and ( fTBI ) together yield 



5 



„ bc + ad , n ac + bd ,„, N 

A 1= 2w — , A 2 = 2 W — — -. (21) 

etc + oa oc + aa 



Note that we kept only positive values of Lagrange multipliers and omitted negative values to get the maximal value of A 
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Now Eq.© takes the form 



max ' 



2 | 8(ab + cd)(ac + bd)(ad + bc)-r 1 r 2 r 3 

In fact, entanglement eigenvalue is the sum of two equal terms and this statement follows from the identity 

r\T2Tz (ab + cd)(ac + bd)(ad + be) 



4uj 2 



(23) 



To derive this identity one has to use the normalization condition a 2 + b 2 + c 2 + d 2 = 1. The identity allows to rewrite Eq.d22l 
as follows 



A max = 4J& (24) 



where 



2 (aft + cd)(ac + bd)(ad + be) 
R « A^r 2 ' (25) 

Above formula has a geometric interpretation and now we demonstrate it. Let us define a quantity p = (a + b + c + d)/2. 
Then the denominator can be rewritten as 

Alj 2 - r| = 16(p - a)(p - b){p - c)(p - d). (26) 

Five independent parameters are necessary to construct a convex quadrangle. However, four independent parameters are 
necessary to construct a convex quadrangle that has circumradius. For such quadrangles the area S q is given exactly by Eq.d26t 
up to numerical factor, that is S 2 = (p — a) (p — b)(p — c) (p — d). Hence Eq.(l25b can be rewritten as 

2 {ah + cd) (ac + bd) (ad + be) 
R « = IRSf ' (2?) 

Thus R q can be interpreted as a circumradius of the convex quadrangle. Eq.(l27b is the generalization of the corresponding 
formula of Ref.[25] and reduces to the circumradius of the triangle if one of parameters is zero. 

Eq. (t24l > is valid if vectors u and v are unit and have non-vanishing x components. These conditions have short formulations 

Kl < 1, \v k \ < 1. (28) 

Above inequalities are polynomials of degree six and algebraic solutions are unlikely. However, it is still possible do define 
the domain of validity of Eq.(T2~7l). 

B. Circumradius of Crossed-Quadrangle 



Here, we consider the second case given by Eq.dTBT). Derivations repeat steps of the previous subsection and the only difference 
is the interchange u> fi. Therefore we skip some obvious steps and present only main results. Components of vectors u and 
v along axis z are 
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u _ A 2 rx - r 2 r 3 ^ _ X x r 2 ~ nr 3 ^ 
— r\ ' 4/i 2 — r\ 

The second equation for Lagrange multipliers 

A x (V 2 + r\ - rl) = A 2 (4 M 2 + r\ - r 2 ) (30) 

together with Eq.(fl~6b yields 

„ be — ad t n ac — bd 

A x =±2/i -, A 2 = ±2/x- -. (31) 

ac — bd be — ad 

Using these expressions, one can derive the following expression for entanglement eigenvalue 

4A 2 =1 A 2 (V+r 2 -r 2 )-r 1 r 2 r 3 

Now the restrictions 1/4 < A 2 lax < 1 derived in Ref.[24] uniquely define the signs in Eq.PTli. Right signs enforce strictly 
positive fraction in right hand side of Eq.d32l. To make a right choice, we replace d by — d in the identity d23l and rewrite Eq.d32b 
as follows 

^ 2 1 (ac — bd){bc — ad){ab — cd) 1 (ac — bd)(bc — ad)(ab — cd) 

max 2 p(p — c — d) (p — b — d) (p — a — d) 2 p(p — c — d)(p — b — d)(p — a — d) 

Lower sign yields zero and is wrong. It shows that reduced density matrix p AB still has zero eigenvalue. 
Upper sign may yield a true answer. Entanglement eigenvalue is 

AL X = 4R 2 X , (34) 

where 

2 (ac - bd)(bc - ad)(ab - cd) 

R * 165 2 ' (35) 

and 5 2 = p(p — c — d)(p — b — d)(p — a — d) . The formula (|35T > may seem suspicious because it is not clear whether right hand 
side is positive and lies in required region. To clarify the situation we present a geometrical treatment of Eq.(f35T>. 

The geometrical figure ABCD in Fig.lA is not a quadrangle and is not a polygon at all. The reason is that it has crossed sides 
AD and BC. We call figure ABCD crossed-quadrangle in a figurative sense as it has four sides and a cross point. Another 
justification of this term is that we will compare figure ABCD in Fig. 1 A with a convex quadrangle ABCD containing the same 
sides. 

Consider a crossed-quadrangle ABCD with sides AB = a, BC = b, CD = c, DA = d that has circumcircle. It is easy to 
find the length of the interval AC 

AC* = (aC ~ M)(6c ~ ad) . (36) 
ab — cd 

This relation is true unless triangles ABC and ADC have the same height and as a consequence equal areas. Note that S x is 
not an area of the crossed-quadrangle. It is the difference between the areas of the noted triangles. 
Using Eq.d36l). one can derive exactly Eq.d35t for the circumradius of the crossed-quadrangle. 
Eq.d34l is meaningful if vectors u and v are unit and have nonzero components along the axis y. 
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Fig. 1A Fig. 1B 



FIG. 1: This figure shows the example for the case when crossed quadrangle(Fig. 1A) has larger circumradius than that of convex quadran- 
gle(Fig.lB) with same sides. 



C. Largest Coefficient 

In this subsection we consider the last case described by Eq.(fT71). Entanglement eigenvalue takes maximal value if all terms 
in r.h.s. of Eq.© are positive. Then equations ( fT71 ) and ( fTOb together impose 

u = Sign(r x )fe, v = Sign(r 2 )fc, nr 2 r 3 < 0, (37) 

where Sign(x) gives -1, or 1 depending on whether x is negative, zero, or positive. Substituting these values into Eq.Q, we 
obtain 



ALx^a + H + M + N). 08) 

Owing to inequality, rir 2 r 3 < 0, above expression always gives a square of the largest coefficient I 

I = max(a, b, c, d) (39) 

in Eq.®. Indeed, let us consider the case r\ > 0,r 2 > 0,r 3 < 0. From inequalities r\ > 0,r 2 > it follows that c 2 > 
d 2 + | a 2 — b 2 \ and therefore c 2 > d 2 . Note, c 2 > d 2 is necessary but not sufficient condition. Now if d > b, then r\ > yields 
c > a and if d < b, then r% < yields c > a. Thus inequality c > a is true in all cases. Similarly c > b and c is the largest 
coefficient. On the other hand = c 2 and Eq.(l38b really gives the largest coefficient in this case. 

Similarly, cases r-y > 0, r 2 < 0, r$ > and r\ < 0, r 2 > 0, > yield A 2 lax = b 2 and A 2 lax = a 2 , respectively. And again 
entanglement eigenvalue takes the value of the largest coefficient. 

The last possibility r\ < 0, r 2 < 0, r$ < can be analyzed using analogous speculations. One obtains A 2 nax = d 2 and d is 
the largest coefficient. 

Combining all cases mentioned earlier, we rewrite Eq.d38l) as follows 



Af nax = F. (40) 
This expression is valid if both vectors u and v are collinear with the axes z. 

We have derived three expressions for (l24ii,(|34"l) and d40b for entanglement eigenvalue. They are valid when vectors u and v 
lie in xz plane, lie in yz plane and are collinear with axis z, respectively. The following section goes on to specify these domains 
by parameters a, b, c, d. 
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V. APPLICABLE DOMAINS 

Mainly, two points are being analyzed. First, we probe into the meaningful geometrical interpretations of quantities R q and 
R x . Second, we separate validity domains of equations d24ii,(|34"ii and d40i i. It is mentioned earlier that algebraic methods for 
solving the inequalities of degree six are ineffective. Hence, we use geometric tools that are elegant and concise in this case. 

We consider four parameters a, b, c, d as free parameters as the normalization condition is irrelevant here. Indeed, one can 
use the state \ip) / %/a 2 + b 2 + c 2 + d 2 where all parameters are free. If one repeats the same steps, the only difference is that 
the entanglement eigenvalue is replaced by A 2 lax / (a 2 + b 2 + c 2 + d 2 ). In other words, normalization condition re-scales 
the quadrangle, convex or crossed, so that the circumradius always lies in the required region. Consequently, in constructing 
quadrangles we can neglect the normalization condition and consider four free parameters a,b,c, d. 

A. Existence of circumcircle. 

It is known that four sides a, b, c, d of the convex quadrangle must obey the inequality p — I > 0. Any set of such parameters 
forms a cyclic quadrilateral. Note that the quadrangle is not unique as the sides can be arranged in different orders. But all these 
quadrangles have the same circumcircle and the circumradius is unique. 

The sides of a crossed-quadrangle must obey the same condition. Indeed, from Fig.lA it follows that BC — AB < AC < 
AD + DC and DC — AD < AC < AB + BC. Therefore AB + AD + DC > BC and AB + BC + AD > DC. The sides 
BC and DC are two largest sides and consequently p — I > 0. However, the existence of the circumcircle requires an additional 
condition and it is explained here. The relation r% — 2/i cos ABC forces 4/x 2 > r 2 and, therefore 

S 2 X > 0. (41) 

Thus the denominator in Eq.(l35ll must be positive. On the other hand the inequality AC 2 > forces a positive numerator of the 
same fraction 

(ac - bd) {be - ad) (ab -cd)>0. (42) 

These two inequalities impose conditions on parameters a,b,c, d. For the future considerations, we need to write explicitly the 
condition imposed by inequality (l42l . The numerator is a symmetric function on parameters a, b, c, d and it suffices to analyze 
only the case a > b > c > d. Obviously (ac — bd) > 0, (ab — cd) > and it remains the constraint be > ad. The last inequality 
states that the product of the largest and smallest coefficients must not exceed the product of remaining coefficients. Denote by 
s the smallest coefficient 

s = min(a, b, c, d). (43) 

We can summarize all cases as follows 

l 2 s 2 < abed. (44) 

This is necessary but not sufficient condition for the existence of R y . The next condition S 2 , > we do not analyze because 
the first condition d44l suffices to separate the validity domains. 

B. Separation of validity domains. 

In this section we define applicable domains of expressions (I24l.d34l and (|40T > step by step. 

a. Circumradius of convex quadrangle. First we separate the validity domains between the convex quadrangle and the 
largest coefficient. In a highly entangled region, where the center of circumcircle lies inside the quadrangle, the circumradius 
is greater than any of sides and yield a correct answer. This situation is changed when the center lies on the largest side of the 
quadrangle and both equations ( f24b and ( f40b give equal answers. Suppose that the side a is the largest one and the center lies on 
the side a. A little geometrical speculation yields 



a 2 = b 2 + c 2 + d 2 + 2 — . 

a 



(45) 
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From this equation we deduce that if a 2 is smaller than r.h.s., i.e. 



a 2 < 6 2 + c 2 +d 2 + 2 ^ ; (46) 

a 

then the circumradius-formula is valid. If a 2 is greater than r.h.s in Eq.d45l>. then the largest coefficient formula is valid. The 
inequality d46l ) also guarantees the existence of the cyclic quadrilateral. Indeed, using the inequality 

bc + cd + bd > 3 — , (47) 
a 

one derives 

(b + c + d) 2 >b 2 + c 2 + d 2 + — >a 2 . (48) 

a 

Above inequality ensures the existence of a convex quadrangle with the given sides. 

To get a confidence, we can solve equation u k — ±1 using the relation d45l ). However, it is more transparent to factorize it as 
following: 

(4c 2 - r 2 )(l + u k ) = ™- (b 2 + c 2 + d 2 + 2 -^-a 2 )(a 2 + b 2 + c 2 + 2 -^- d 2 ) (49a) 
bc + ad \ a J \ d J 

(4c 2 - r 2 )(l - u k ) = (a 2 + c 2 + d 2 + 2 -^- b 2 ) ( a 2 + b 2 + d 2 + ™* - c 2 ) . (49b) 



be + ad 



Similarly, we have 



(4c 2 - r 2 )(l + v k ) = (a 2 + c 2 + d 2 + ^ - b 2 ) (a 2 + b 2 +c 2 + ^ -d 2 ) (50a) 

ac + bd \ b I \ d 1 



(4c 2 - r 2 )(l - v k ) = (b 2 + c 2 + d 2 + — - a 2 ) ( a 2 + b 2 + d 2 + — - c 2 ) . (50b) 

ac + bd \ a J \ c J 

Thus, the circumradius of the convex quadrangle gives a correct answer if all brackets in the above equations are positive. In 
general, Eq.(l24l> is valid if 

, 9 1 abed 

l 2 <~ 2 + — (5D 



When one of parameters vanishes, i.e. abed = 0, inequality dBTl ) coincides with the corresponding condition in Re f. 02511 . 

b. Circumradius of crossed quadrangle. Next we separate the validity domains between the convex and the crossed quad- 
rangles. If S* 2 < 0, then crossed one has no circumcircle and the only choice is the circumradius of the convex quadrangle. If 
S 2 , > 0, then we use the equality 



o o r abed 
4i? 2 -4i? 2 =-— -y (52) 



where r = r^r^. It shows that r > yields R q > i? x and vice-versa. Entanglement eigenvalue always takes the maximal 
value. Therefore, = 4i? 2 if r > and A 2 nax = 4i? 2 if r < 0. Thus r = is the separating surface and it is necessary to 

analyze the condition r < 0. 

Suppose a > b > c > d. Then and r% are positive. Therefore r is negative if and only if r\ is negative, which implies 



a 2 + d 2 > b 2 + c 2 . 



(53) 
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Now suppose a > d > b > c. Then n is negative and r% is positive. Therefore r<i must be positive, which implies 

a 2 + c 2 >b 2 + d 2 . (54) 

It is easy to see that in both cases left hand sides contain the largest and smallest coefficients. This result can be generalized 
as follows: r < if and only if 

I 2 > - - s 2 . (55) 
- 2 

It remains to separate the validity domains between the crossed-quadrangle and the largest coefficient. We can use three 
equivalent ways to make this separation: 

1) to use the geometric picture and to see when 4i? x and I 2 coincide, 

2) directly factorize equation Uk = ±1, 

3) change the sign of the parameter d. 

All of these give the same result stating that Eq.(l34b is valid if 

2 1 abed 

I <~ 2 -—- (56) 

Inequalities d55b and d56b together yield 

l 2 s 2 > abed. (57) 

This inequality is contradicted by d44b unless l 2 s 2 — abed. Special cases like l 2 s 2 = abed are considered in the next section. 
Now we would like to comment the fact that crossed quadrangle survives only in exceptional cases. Actually crossed case can 
be obtained from the convex cases by changing the sign of any parameter. It crucially depends on signs of parameters or, in 
general, on phases of parameters. On the other hand all phases in Eq.® can be eliminated by LU-transformations. For example, 
the phase of d can be eliminated by redefinition of the phase of the state function and the phases of remaining parameters 
can be absorbed in the definitions of basis vectors |1) of the qubits A, B and C. Owing to this entanglement eigenvalue being 
LU invariant quantity does not depend on phases. However, crossed case is relevant if one considers states given by Generalized 
Schmidt Decomposition(GSD) BlOll . In this case phases can not be gauged away and crossed case has its own range of definition. 
This range has shrunk to the separating surface r = in our case. 

Now we are ready to present a distinct separation of the validity domains: 



a I ARl if I 2 < 1/2 + abed! 'I 2 
max l I 2 if I 2 > 1/2 + abcd/l 2 



As an illustration we present the plot of ^-dependence of A„ lax in Fig. 2 when a = b = c. 

We have distinguished three types of quantum states depending on which expression takes entanglement eigenvalue. Also 
there are states that lie on surfaces separating different applicable domains. They are shared by two types of quantum states and 
may have interesting features. We will call those shared states. Such shared states are considered in the next section. 

VI. SHARED STATES. 

Consider quantum states for which both convex and crossed quadrangles yield the same entanglement eigenvalue. Eq.(l36i> is 
not applicable and we rewrite equations (|27| > and (l35l > as follows 

^K 1- !^)' 4 ^=K 1- ^)' (59) 

These equations show that if the state lies on the separating surface r = 0, then entanglement eigenvalue is a constant 
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(4/7)" 



FIG. 2: Plot of d-dependence of A^ aa; when a — b — c. When d — * 1, A„ lax goes to 1 as expected. When d = 0, A„ lax becomes 4/9, which 
coincides with the result of Ref.[7]. When r = which implies a = d — 1/2, A^ ax becomes 1/2 (it is shown as dotted line). When d = 2a, 
which implies d = a/4/7, A^ ax goes to 4/7, which is one of shared states (it is also shown as another dotted line). 



A 2 = - 

max 2 



(60) 



and does not depend on the state parameters. This fact has a simple interpretation. Consider the case r\ = 0. Then b 2 + c 2 = 
a 2 + d 2 = 1/2 and the quadrangle consists of two right triangles. These two triangles have a common hypotenuse and legs b, c 
and a, d, respectively, regardless of the triangles being in the same semicircle or in opposite semicircles. In both cases they yield 
same circumradius. Decisive factor is that the center of the circumcircle lies on the diagonal. Thus the perimeter and diagonals 
of the quadrangle divide ranges of definition of the convex quadrangle. When the center of circumcircle passes the perimeter, 
entanglement eigenvalue changes-over from convex circumradius to the largest coefficient. And if the center lies on the diagonal, 
convex and crossed circumradiuses become equal. 

We would like to bring plausible arguments that this picture is incomplete and there is a region that has been shrunk to the 
point. Consider three-qubit state given by GSD 



\ip) = a\W0) + 6|010) + &|001> + + e|000). 



(61) 



One of parameters must have non-vanishing phase[10I] and we can treat this phase as an angle. Then, we have five sides 
and an angle. This set defines a sexangle that has circumcircle. One can guess that in a highly entangled region entanglement 
eigenvalue is the circumradius of the sexangle. However, there is a crucial difference. Any convex sexangle contains a star type 
area and the sides of this area are the diagonals of the sexangle. The perimeter of the star separates the convex and the crossed 
cases. Unfortunately, we can not see this picture in our case because the diagonals of a quadrangle confine a single point. It is 
left for future to calculate the entanglement eigenvalues for arbitrary three qubit states and justify this general picture. 

Shared states given by r = acquire new properties. They can be used for perfect teleportation and superdense coding 
l25l I27I l28ll . This statement is not proven clearly, but also no exceptions are known. 

Now consider a case where the largest coefficient and circumradius of the convex quadrangle coincide with each other. The 
separating surface is given by 



abed 



(62) 



Entanglement eigenvalue ranges within the narrow interval 



12 




max — 



4 

< - 



- 7' 



(63) 



It separates slightly and highly entangled states. When one of coefficients is large enough and satisfies the relation I 2 > 
1/2 + abcd/l 2 , entanglement eigenvalue takes a larger coefficient. And the expression ([8]) for the state function effectively takes 
the place of Schmidt decomposition. In highly entangled region no similar picture exists and all coefficients participate in equal 
parts and yield the circumradius. Thus, shared states given by Eq.d62l) separate slightly entangled states from highly entangled 
ones, and can be ascribed to both types. 

What is the meaning of these states? Shared states given by r = acquire new and important features. One can expect that 
shared states dividing highly and slightly entangled states also must acquire some new features. However, these features are yet 
to be discovered. 



We have considered three-parametric families of three qubit states and derived explicit expressions for entanglement eigen- 
value. The final expressions have their own geometrical interpretation. The result in this paper with the results of Ref.[25] show 
that the geometric measure has two visiting cards: the circumradius and the largest coefficient. The geometric interpretation 
may enable us to predict the answer for the states given by GSD. If the center of circumcircle lies in star type area confined 
by diagonals of the sexangle, then entanglement eigenvalue is the circumradius of the crossed sexangle(s). If the center lies in 
the remaining part of sexangle, the entanglement eigenvalue is the circumradius of the convex sexangle. And when the center 
passes the perimeter, then entanglement eigenvalue is the largest coefficient. Although we cannot justify our prediction due to 
lack of computational technique, this picture surely enables us to take a step toward a deeper understanding of the entanglement 
measure B29I1 . 

Shared states given by r — play an important role in quantum information theory. The application of shared states given 
by Eq.d62b is somewhat questionable, and should be analyzed further. It should be pointed out that one has to understand the 
properties of these states and find the possible applications. We would like to investigate this issue elsewhere. 

Finally following our procedure, one can obtain the nearest product state of a given three-parametric W-type state. These 
two states will always be separated by a line of densities composed of the convex combination of W-type states and the nearest 
product states [30]. There is a separable density matrix qq which splits the line into two parts as follows. One part consists 
of separable densities and another part consists of non-separable densities. It was shown in Ref. 113011 that an operator W — 
go — p ABC — tr[go(go — p ABC )]I has the properties tr(Wp ABC ) < 0, and tr(Wg) > for the arbitrary separable state g. The 
operator W is clearly Hermitian and thus is an entanglement witness for the state. Thus our results allow oneself to construct 
the entanglement witnesses for W-type three qubit states. However, the explicit derivation of go seems to be highly non-trivial 
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